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The proposal of a galactic population of intermediate mass black holes (IMBHs), forming dark 
matter (DM) "mini-spikes" around them, has received considerable attention in recent years. In 
fact, leading in some scenarios to large annihilation fluxes in gamma rays, neutrinos and charged 
cosmic rays, these objects are sometimes quoted as one of the most promising targets for indirect 
DM searches. In this letter, we apply a detailed statistical analysis to point out that the exist- 
ing EGRET data already place very stringent limits on those scenarios, making it rather unlikely 
that any of these objects will be observed with, e.g., the Fermi/GLAST satellite or upcoming Air 
Cherenkov telescopes. We also demonstrate that prospects for observing signals in neutrinos or 
charged cosmic rays seem even worse. Finally, we address the question of whether the excess in the 
cosmic ray positron/electron flux recently reported by PAMELA/ ATIC could be due to a nearby 
DM point source like a DM clump or mini-spike; gamma-ray bounds, as well as the recently released 
Fermi cosmic ray electron and positron data, again exclude such a possibility for conventional DM 
candidates, and strongly constrain it for DM purely annihilating into light leptons. 
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The nature of the mysterious DM, vastly dominating 
the total matter content of the universe, still remains un- 
known. Particularly plausible candidates, however, are 
weakly interacting massive particles (WIMPs) and in- 
direct DM searches aim at discriminating WIMP annihi- 
lation products from standard astrophysical backgrounds 
in gamma rays, neutrinos or charged cosmic rays. 

Large DM density enhancements ("mini-spikes") 
around IMBHs have been proposed as promising targets 
for indirect DM searches in gamma rays , where a large 
number of very luminous point sources with identical cut- 
off in the photon spectrum would provide a smoking gun 
signature. Subsequent studies indicated excellent obser- 
vational prospects also for neutrinos Q and charged cos- 
mic rays In the following, we restrict ourselves to 
IMBHs with a mass of around lO^M© that form out of 
collapsing cold gas in early-forming halos (scenario B of 
Ref. i), Hke in much of the literature on the subject, 
since only in this case one arrives at the above mentioned 
favorable prospects for indirect DM detection that have 
caused considerable recent attention. 

While an application of our analysis to other IMBH for- 
mation scenarios would be straightforward, it is beyond 
the scope of this Letter, in which we critically re-assess 
the potential of IMBHs for DM searches. Our main con- 
clusion is that the most favored DM parameter regions 
are actually already ruled out by the EGRET data [5| and 
that configurations predicting a signal in future searches 
- and yet being consistent with the existing constraints ~ 



are rather unlikely. For completeness, we will also derive 
generic limits on any nearby DM point-source. 

Let us start by recalling some basic statistical proper- 
ties of an ensemble of IMBHs in the galactic halo. De- 
noting with Pc the probability for a particular IMBH to 
satisfy some condition c, the probability for n out of N 
objects to satisfy c is given by 



Pn,N{Pc) = 



P^ (1 - Pc) 



N-r, 



(1) 



Since the number of IMBHs in a given realization itself 
is a random variable with some distribution p^v (in our 
case a Gaussian with mean ~ 98 and variance ~ 21 0]), 
the probability that n objects in an arbitrary realization 
satisfy c becomes 



dNpN{N)P„MiPc)- 



(2) 



In our context, the relevant astrophysical properties 
of an IMBH are its distance d to the Earth and, as a 
measure of the DM concentration, its annihilation volume 
£, = J <i^x (p(x)/po)^, where po = O.SGeVcm"'^ is the 
local DM density. In particular, we will be interested in 
the probability that a single object has A = ^/d^ > Ao: 

Pa>a„ = / / dd{p^{OPd{d)), (3) 

Jo Jo 

where we use the (independent) probability densities pd 
and p^ to find an IMBH at a distance d from the earth 



2 



and with a given annihilation volume ^, respectively, from 
Ref. [1]. These distributions were obtained from Monte 
Carlo (MC) simulations that follow the evolution of ini- 
tial mini-spike populations during their orbit in the Milky 
Way, taking account of possible close encounters. Armed 
with the above notation, we finally arrive at 
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n-l 

E 

1=0 



(4) 



as the probability that a given realization contains at 
least n objects with A > Aq. 

The gamma-ray flux from a single IMBH is given by 



^IMBH 



(5) 



where F = ■^{(Jv) (po/"^x)^^ ^'^^ local injection rate, 
771^ the WIMP mass, (av) the annihilation rate (per unit 
density) and the number of photons above some en- 
ergy -Eniin per annihilation. At <C m^, the pho- 
ton spectrum looks almost exactly the same for generic 
WIMP DM candidates that annihilate into quark or 
gauge boson final states; if not stated otherwise, we as- 
sume annihilation into bb and use Pythia Q to compute 
N-y. The adiabatic growth of an IMBH redistributes a 
typical initial DM density distribution into a steep profile 
(a " mini-spike" ) that is saturated in the innermost region 
due to DM self-annihilations and asymptotically develops 
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an annihilation volume that scales like ^ oc (m^/ {av)) ' 
. Putting in numbers, Eq.([5]) thus becomes 



i 



AT 

^IMBH 



3.31 




(6) 



(av) 



3-10- 



where ^ is the annihilation volume for an IMBH when 
{av) = 3 • 10-2^cm3/s and = 1 TeV. 

The upper limit on the flux from point sources not 
seen by the EGRET satellite strongly depends on the 
position on the sky, but is almost everywhere significantly 
below - 2 X lO-'^cm-^s-i for 100 MeV < < 

30GeV [51 . On the other hand, a considerable number 
of detected sources with higher fluxes remains with no 
associated low-energy counterpart. Very conservatively 
assuming that all these unidentified sources arc, in fact, 
connected to DM mini-spikes would translate into the 
requirement that at most (1, 4, 10) IMBHs have a flux 
larger than about $L^ax^°^ ~ (7, 5.5, 3.7) x 10"^ cm'^s"! 
; as the spectra of these sources are in almost all cases 
much softer than what is expected from DM annihilation, 
however, this possibility does not appear to be very likely 
(note also that many of these EGRET sources are not 
confirmed by the first Fermi data Q). 

Demanding that at least 5% of the IMBH realizations 
should not be in conflict with the EGRET constraints 
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FIG. 1: The thin solid lines show the constraints on WIMPs 
in the IMBH scenario B of 0], from top to bottom deriving 
from the brightest 10, 4 and 1 object(s), respectively. The 
region above the thick solid line is excluded by the EGRET 
limit on unresolved point sources; the dotted lines show the 
corresponding constraint if DM were to annihilate only into 
e^e~ or fj,'^ . The dashed line indicates the canonical value 
of av = 3 • lO"^'^ cm^s"^ for thermally generated WIMPs. 



on the n brightest sources, we can now use Eq. @ to 
derive bounds on the annihilation rate and mass of the 
annihilating particles by choosing Xq = ^/(P such that 
P™™(Ao) > 0.95. The result is shown in Fig. [T] (for simi- 
lar constraints from the H.E.S.S. experiment, see Q). It 
is interesting to note that the n ^ 10 brightest unidenti- 
fied EGRET sources lead to similar constraints. Having 
remarked before, however, that an IMBH interpretation 
of these objects is not too likely, we also include for com- 
parison the much tighter constraint that results from con- 
sidering the sensitivity limit of EGRET on unseen point 
sources. We verified these limits in extensive MC sim- 
ulations and would like to draw special attention to the 
fact that they lie, indeed, several orders of magnitude 
below the expectation for generic WIMP candidates [l[. 



0(100 GcV-1 TeV) and 



0(10 



cm^s-i). 



Note that the constraints shown in Fig. [T] are, in fact 
even conservative as the EGRET limit on unresolved 
point sources is actually in many regions of the sky much 
smaller than the value used here. 

The neutrino fiux is simply obtained by replacing 
in Eq. (O the index 7 with ly, so the EGRET hmit 
on gamma rays translates into 'I'Jmbh < (^V-^'')*J'max- 
This limit on the neutrino flux has to be compared with 
the sensitivity of upcoming km'^-sized neutrino telescopes 
[Tot . Optimistically assuming an effective surface area 
for detection below 1 TeV of - 0.l{E^/TeV)^m'^, and us- 
ing the Bartol model for the atmospheric neutrino back- 
ground Uj, we find a background rate of ~ 10~^ Hz. 
A 5(7 detection after 1 yr would then require a primary 
neutrino fiux 'I'imbh ^ 10~®cm~2s~^ (over-optimistically 
assuming that the target is observed 100% of the time). 
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which exceeds the above stated hmit for typical values 
of rrij^ and N'^. Of course, should DM annihilate mostly 
into neutrinos, the EGRET limit would be less stringent. 
However, even in conventional Kaluza-Klein (KK) sce- 
narios [12| . which offer large branching ratios to neutri- 
nos, the EGRET constraint still excludes the observation 
of the neutrino counterpart with km'^ detectors. Note 
also that large branching ratios into neutrinos actually 
mean smaller TV^, at least for large m^, leading to over- 
all worse prospects for detection [3|]. 

Complementary to gamma rays, antiprotons are a fur- 
ther interesting channel of indirect DM detection . To 
investigate whether the resulting constraints can compete 
with those from gamma rays in the mini-spike scenario, 
we have performed an extensive MC simulation of IMBH 
realizations. We found that antiprotons are competitive 
(i.e. in potential conflict with high-energy p data 14{ ) 



only for rather heavy WIMPs, with rii^ > 1 TeV, and 
only when assuming extremely favorable pro pag ation pa- 
rameters, close to the max set defined in [l5[; outside 
these somewhat extreme regions of the propagation pa- 
rameter space, large p fluxes are generally obtained only 
in configurations that are already excluded by EGRET. 

The excess in cosmic ray positrons and electrons re- 
cently reported by PAMELA ATIC [Hi and Fermi 
[ist . if interpreted in terms of DM annihilation, points 
at DM masses in the TeV range with an unusually large 
branching ratio into light leptons • In order to fit the 
data, however, standard DM candidates need extremely 
large boost factors that are not expected in current mod- 
els of structure formation [2^ , corresponding to effective 
annihilation rates 10^ times the generic value for ther- 
mally produced DM. For comparison, we thus further- 
more include in Fig. [T] the extreme situation of DM par- 
ticles annihilating only into light leptons, in which case 
much less photons are produced and we can apply the an- 
alytic expression for given in 2]|; even in this some- 
what contrived situation low-mass models turn out to be 
difficult to realize, at least for strongly enhanced annihi- 
lation rates, making it very unlikely for, e.g, Fermi-LAT 
[2^ to see the expected cutoff in the photon spectrum 
for an IMBH that has escaped detection by EGRET. 

Let us now turn to the possibility that the large re- 
quired annihilation flux can be attributed to a nearby 
high DM concentration, like for example a DM clump 
or a DM mini-spike around an IMBH, located at some 
distance d to the Earth. In Fig. [21 we use Eq. ([5]) and 
the EGRET sensitivity limit to constrain the annihila- 
tion rate F of such a generic DM point source. In the 
same figure, we indicate as a gray area in the T — d 
plane the annihilation rate that would be necessary to ex- 
plain the PAMELA data, taking into account the allowed 
range for the e^/e~ propagation parameters (consistent 
with both the background and the signal) by using the 
min/med/max configurations of psj . The dark shaded 
region, finally, shows the rather conservative constraint 



on these combinations of F and d that arises from re- 
quiring that the positron and electron flux from DM an- 
nihilation alone should not exceed the Fermi data. For 
comparison, we separately consider the case of KK DM, 
as an example of a standard WIMP with exceptionally 
large branching ratios into light leptons, as well as DM 
only annihilating to e~^e~ or 

The first important observation is that standard 
WIMPs, which usually feature smaller branching ratios 
into light leptons than KK DM, basically cannot account 
for the PAMELA/ ATIC data in this way without vio- 
lating the EGRET bounds - the reason being the still 
relatively large contribution from non-leptonic channels 
to the photon spectrum at low energies. However, if one 
takes a more phenomenological approach and allows DM 
particles annihilating at 100% into e^e~ or /i^^/i" pairs, 
one would technically be able to fit the data by plac- 
ing a DM point source at a distance 2 kpc ^ d < b kpc 
- at the price of requiring an enormously bright object, 
more luminous than the whole Milky Way\ This would 
correspond to a DM clump of mass > IO^^Mq in con- 
ventional cosmological scenarios [2^ . Since finding such 
a massive clump relatively close to the Earth is extremely 
unlikely , we arrive at a considerably more pessimistic 
conclusion than recently obtained by [25| about the pos- 
sibility of explaining the PAMELA data in terms of a 
nearby clump of annihilating WIMPs. Mini-spikes, on 
the other hand, are extremely bright: for DM with 
{(jv) ^ 3 • lO^^^cm'^s"^ in the IMBH scenario discussed 
earlier, e.g., the probability to encounter at least one ob- 
ject inside the light gray area of the middle (right) panel 
of Fig. [2 is roughly 84% (37%). For a - certainly non- 
standard - WIMP candidate annihilating almost exclu- 
sively into e+e~ or /x^/i", a galactic population of IMBHs 
might thus indeed provide a positron flux large enough 
to fit the data and yet be consistent with present-day 
constraints. Note, however, that this conclusion does not 
hold for DM candidates with intrinsically enhanced an- 
nihilation rates: already for {av) ~ lO^^^cm'^s^^, the 
above quoted probabihties drop to 0.9% (0.1%). 

In this Letter, we have reconsidered the prospects for 
indirect DM detection in the presence of a galactic popu- 
lation of IMBHs and found them not to be very promising 
given that existing data, in particular from gamma rays, 
already place severe constraints on the scenario. While 
it was noted before that already EGRET should have 
seen some of these objects 0, the resulting constraints 
were not taken into account in the subsequent studies of, 
e.g., [3, 0]- This is the main reason for the discrepancy 
between our pessimistic and earlier rather optimistic con- 
clusions about possible effects of DM mini-spikes on in- 
direct DM searches. While beyond the scope of this let- 
ter, it would be interesting to extend the study presented 
here and apply a consistent treatment of all available con- 
straints also to other IMBH formation scenarios in order 
to predict realistic prospects for indirect DM detection. 
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FIG. 2; The solid lines give the EGRET constraints on the DM annihilation rate T = ^crt; (po/wix)^ C of ^ nearby, generic 
DM point-source at a distance d from the Earth; from left to right, we show the case of KK DM and a fiducial DM candidate 
annihilating to e~ and n'^ fj,~ , respectively. The dashed lines show the F needed to fit the PAMELA data, for sets of 
propagation parameters as defined in [SJ; in the dark shaded area this would produce an flux in conflict with the Fermi 
data at higher energies. For comparison, the dotted line indicates F for the whole Milky Way, assuming (av) ~ 3-10~'^^cm"''s^^. 



We have also addressed the possibility of explaining 
the recent PAMELA observations in terms of DM anni- 
hilation by placing a dark object in close vicinity to the 
Earth. While this option is ruled out from the gamma- 
ray constraints for standard WIMP candidates, a nearby 
IMBH (but not an ordinary DM clump) may in principle 
provide a sufficient amount of positrons if one assumes 
that the DM particles annihilate purely into light Icptons. 
Even this seemingly a bit far-fetched scenario could soon 
be ruled out if Fermi [2^ does not observe any corre- 
sponding point-sources in gamma rays. We would like 
to use this opportunity to recall that a very plausible 
hypothesis for the PAMELA results is anyway a nearby 
pulsar, i.e. an explanation in terms of astrophysics rather 
than DM annihilation [2^. 

Acknowledgements: We greatly thank Gianfranco 
Bertone for letting us use his catalogue of IMBH sources 
and the anonymous referees for very useful suggestions. 
TB and JL thank LAPTH for hospitality when this work 
was initiated. This work is supported, in part, by the 
French ANR project ToolsDMCoU (BLAN07-2-194882). 
PS thanks the Perimeter Institute for hospitality when 
this work was finished. 



* Electronic address: troms@physto.se 
^ Electronic address: lavalle@to. infn.it 

* Electronic address: salati@la pp.in2p3.fr| 

[1] G. Jungman, M. Kamionkowski and K. Griest, Phys. 

Rept. 267, 195 (1996); L. Bergstrom, Rept. Prog. Phys. 

63, 793 (2000); G. Bertone, D. Hooper and J. Silk, Phys. 

Rept. 405, 279 (2005). 
[2] G. Bertone, A. R. Zentner and J. Silk, Phys. Rev. D 72, 

103517 (2005). M. Fornasa and G. Bertone, Int. J. Mod. 



Phys. D 17, 1125 (2008). 
[3] G. Bertone, Phys. Rev. D 73, 103519 (2006). 
[4] P. Brun et al, Phys. Rev. D 76 (2007) 083506. 
[5] Hartman et al., Ap.J.S. 123, 79 (1999). 
[6] T. Sjost rand, S. Mrenna and P . Skands, JHEP 0605, 026 

(2006) |: arXiv:hep-ph/0603175| . 

[7] V. Pavlidou et aL, Astrophys. Space Sci. 309 (2007) 81. 

[8] A. A. Abdo et al. l arXiv:0902T 340 [astro-ph.HE]. 

[9] F. Aharonian et al, Phys. Rev. D 78, 072008 (2008). 
[10] See, e.g., T. Montaruh, arXiv:0901,26_61 [astro-ph]. 
[11] G. Barr, T. K. Gaisser and T. Stanev, Phys. Rev. D 39 
(1989). 

[12] D. Hooper and S. Profumo, Phys. Rept. 453, 29 (2007). 
[13] T. Bringmann and P. Salati, Phys. Rev. D 75, 083006 

(2007) ; F. Donato et al., Phys. Rev. Lett. 102, 071301 
(2009). 

[14] O. Adriani et al, Phys. Rev. Lett. 102, 051101 (2009). 

[15] F. Donato et al., Phys. Rev. D 69, 06351 (2004). 

[16] O. Adriani et a/., larXiv:0810.4995 [astro-ph]. 

[17] J. Chang et al. Nature 456, 362 (2008). 

[18] A. A. Abdo et al. Phys. Rev. Lett. 102, 181101 (2009). 

[19] see, e.g., M. Cirelh, M. Kadastik, M. Raidal and A. Stru- 
mia, Nucl. Phys. B 813, 1 (2009). 

[20] V. Springel et al, MNRAS 391, 1685 (2008). 

[21] L. Bergstrom, T. Bringmann, M. Eriksson and 
M. Gustafsson, Phys. Rev. Lett. 94, 131301 (2005). 

[22] E. A. Bahz et al, JCAP 0807, 013 (2008). 

[23] T. Delahaye et a/., larXiv:0 809.5268 [astro-ph]. 

[24] J. Lavalle, Q. Yuan, D. Maurin and X. J. Bi, Astron. 
Astrophys. 479, 427 (2008). 

[25] D. Hooper, A. Stebbins and K. M. Zurek, 
larXiv:0812.3202l [astro-ph] . 

[26] A. Boulares, ApJ 342, 807 (1989); F. A. Aharonian, 
A. M. Atoyan and H. J. Voelk, Astron. Astrophys. 294, 
L41-L44 (1995); H. Yiiksel, M. D. Kistler and T. Stanev, 
arXiv:0810.2784; D. Hooper, P. Blasi and P. D. Ser- 
pico, JCAP 0901, 025 (2009); S. Profumo, N. J. Shaviv, 
E. Nakar and T. Piran, iarXiv:0902.0376, [astro-ph.HE]. 



